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Abstract 

Shadowing corrections to the structure functions of heavy nuclei are calcu- 
lated at very low values of Bjorken-x and at values of the momentum transfer 
relevant to recent experiments. Good agreement is obtained with data from 
the E665 Collaboration for Xe/D and Pb/D, and with the NMC data on 
Ca/D and C /D structure function ratios. Corrections to the deuteron struc- 
ture function are also estimated down to x ~ 10~ 5 , and found to be less than 

about 3% over the range of x covered by the E665 data. 
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Recent experiments by the E665 Collaboration at Fermilab have measured deep 
inelastic nuclear structure functions at values of Bjorken-x some two orders of magnitude 
smaller than those attainable in previous experiments. The data on heavy nuclei (Xe and 
Pb) complement the existing high-precision NMC data on lighter nuclei, such as Ca and 
C, which will help us to understand the origin of the nuclear EMC effect at small x. In 
addition, the E665 Collaboration measured the deuteron to proton structure function ratio 
101, from which the neutron structure function was determined down to x ~ 10~ 5 . To date, 
however, the shadowing corrections to i<2 n have not been estimated in the region of x and 
Q 2 covered by the E665 data. We do so here. Furthermore, we calculate the shadowing 
effect in heavier nuclei, which will also serve as a check of the reliability of the deuteron 
predictions. 

We consider a two-phase model which incorporates hadron and parton degrees of freedom 
at low and high Q 2 values, respectively. The approach is to use the vector meson dominance 
(VMD) model, which describes low-Q 2 , virtual-photon interactions, and to model the ap- 
proximate scaling of shadowing by Pomeron-exchange. A synthesis of these two mechanisms 
is quite sensible if we are to investigate the transition region between small- and large-Q 2 
processes. The physical picture is that the virtual photon interacts with the nuclear target 
via its fluctuations into qq pairs. If the virtuality of the photon is large the fluctuation is 
short-lived, and a description in terms of diffractive scattering from partons is appropriate. 
If the virtuality is smaller, the virtual qq pair will have time to evolve into a state resembling 
a vector meson, which then enables a VMD-based description to be used. In this sense the 
approach is similar to the model of shadowing adopted by Kwiecinski & Badelek || , and to 
the recent description of real photon interactions by Schuler & Sjostrand j4|. 

The empirical basis for the VMD model is the observation that photon — hadron pro- 
cesses have many remarkable similarities with purely hadronic reactions. The most simple 
and natural explanation of this phenomenon is that the photon itself has a hadronic (vector 
meson) structure. Since VMD gives a very good description of shadowing in high energy pho- 
toproduction, this suggests that the same mechanism may also be responsible for shadowing 
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in deep inelastic scattering (DIS), at least at low Q 2 . Formally, the amount of shadow- 
ing can be quantified via the Glauber multiple scattering expansion, which, in the eikonal 
approximation, gives a correction to the nuclear structure function (per nucleon): 

A Q>) = ? E <^ (la) 



where 



2 iz'>2 

is the shadowing correction to the V A cross section, with b the impact parameter, k\ = 
MV(1 + M^/Q 2 ) 2 , and M the nucleon mass. By the Heisenberg uncertainty principle 
only the lowest mass vector mesons (p ,u,<f>) will be important, and these couple to the 
photon with strengths fy = 2.28, 26.14 and 14.91, respectively ||. The approximation in 
Eqs.(l) also omits non-diagonal vector meson transitions (VN —± V'N), however these are 
not expected to be large. 

In Eq. (|TT5|) p^(r, r') = N c p(r) p(r') {1 — C(r — r')} is the two-body density function 
(with p(r) the single body density) normalised so that J d 3 r d 3 r' p^(r, r') = / d 3 r p(r) = 1. 
The correlation function C(r — r') takes into account the short-range repulsion of the NN 
force, and can be modelled by 

C(r-rQ=f 3j ' l( , l6|r "f l) V (2) 



K\r 



where k = 3.6 fm _1 is chosen to reproduce a 'hole' in the two-body density which is ~ 0.5 fm 
wide at 1/2 maximum density. For the single particle density in heavy nuclei (A ^ 16) we use 
the Woods-Saxon (or Fermi) density, while for light nuclei (A < 16) an harmonic oscillator 
(shell model) form is more appropriate 0. 

For light nuclei, the dominant process is that involving the double scattering of vector 
mesons from two nucleons. Higher order terms (multiple rescattering) in the Glauber ex- 
pansion attenuate the incident flux of vector mesons as they traverse the nucleus, which will 



be progressively more important as A increases. The magnitude of the attenuation is de- 
termined by the mean free path of the vector meson in the nucleus, Ly = (<Jvn p(b, C)) _1 
[/J. For the total VN cross sections we use the energy dependent forms from Ref. ||: 
(TpOjv = (Jujn = 13.63 s e + 31.79 s~ v and a^N = 10.01 s t — 1.51 s~ v , where s = (p+ q) 2 , with 
p and q the nucleon and photon four-momenta, respectively. The parameters e ~ 0.08 and 
7] « 0.45 are taken from Regge theory. 

Numerically, there is a strong Q 2 dependence of 8^F 2 a, which peaks at around Q 2 ~ 1 
GeV 2 . For Q 2 — > 0, 5^i*2A disappears due to the vanishing of the total F 2 a- Furthermore, 
since this is a higher twist effect (because of the vector meson propagators in Eq . flTa]) ), 
shadowing in the VMD model dies off quite rapidly between Q 2 ~ 1 and 10 GeV 2 , so that 
for Q 2 > 10 GeV 2 it is almost negligible. In order to reproduce the observed persistence of 
shadowing in this region j9], one could extend the model by including additional hadronic 



states, or a high-mass qq continuum, as in generalised vector meson dominance models ||T0 
However, when discussing DIS phenomena at high Q 2 it may be more efficient to use a 
partonic description. This is certainly advantageous when describing scaling of inelastic 
nucleon structure functions in the parton model. 

The parton model description of diffractive deep inelastic processes in the Bjorken limit 
corresponds to the Regge region (s m Q 2 /x 3> Q 2 )- Therefore nucleon DIS at small x can be 
viewed in terms of virtual-photon interactions with the Pomeron (V) through its 'structure 
function', F 2 v- In calculations of F 2 v two contributions are usually included, an explicit 
qq component [p] , p^| , and the triple Pomeron term |l2|fl, F 2 -p(x-p,Q 2 ) = F 2 p\x-p,Q 2 ) + 
F^ixv^Q 2 ), where 0JT3 



fgW, ?) = 4(10 IrfUa-^) (^) 

F^\x v ,Qi) = J!^FS(x v ,Q*). (3b) 

In Eqs.(0) x-p = x/y is the fraction of the Pomeron's momentum carried by the struck 
quark, y(= k ■ q/p ■ q = x(l + M^/Q 2 ) ~ M^/s) is the fraction of the momentum of the 
nucleon carried by the Pomeron (four-momentum k) and M\ = (k + q) 2 is the mass squared 



of the produced hadronic debris X. In Eq.(|3a]) (3q = 3.4 GeV 2 is the quark — Pomeron 



coupling constant [P|| , and X s ~ 0.5 represents the weaker coupling of V to the strange 
quark compared with the u and d quarks. The parameter N sea 0.17 at Q 2 ~ 4 GeV 2 ) 
is determined by the small- x behaviour of the nucleon sea distribution, xq sea (x —>■ 0) — > 
iV^x 1 "^ ), where a P (0) « 1 is the P intercept |TTJ. In Eq.(§b]) g> 37 > = 0.364 mb 1/2 is 



the triple Pomeron coupling constant |14J. For the pp total cross section we use the Regge 
theory-based parameterisation: a pp = 21.70s e + 56.08s _,? ||. 

For fixed z/, the small- X region also corresponds to small Q 2 . At current energies, this 
means that probing nuclear structure functions at x ^ 10~ 3 is possible only for Q 2 < 0.5 
GeV 2 , which is well below the scaling region. Since the Q 2 dependence of F 2 -p is determined 
by the nucleon sea quark densities, we need to know how these behave at low Q 2 . Although a 
rigorous theoretical basis for nucleon structure functions at very small Q 2 is still outstanding, 
some phenomenological parameterisations have been constructed [f4|,|15"! by incorporating 
the photoproduction limit. The simplest, and most common method adopted has been to 
multiply F*2N by a factor (Q 2 /(Q 2 + Qq)) 1+6 , with the normalisation and mass parameter 
Ql fixed by matching the photoproduction and deep inelastic regions. In Eq. ( ]3lo| ) we extract 
from the parameterisation of F 2p and F 2n in Ref. and for the qq contribution we 



include in Eq.([3a|) the factor Q 2 /(Q 2 + Ql)- The parameter Ql « 0.485 GeV 2 , which we take 
from Donnachie & Landshoff |15| , may be interpreted as the inverse size of partons inside 
the virtual photon. 

In diffractive DIS from nuclei Pomeron exchange between the projectile and two or more 
constituent nucleons models the interaction of partons from different nucleons within the 
nucleus. The shadowing correction to F 2 a from "P-exchange can be written as a convolution 
of F 2 -p and an exchanged-P distribution function, f-p/A- 

A 5^F 2A (x, Q 2 ) = [ A dy f v/A (y) F 2P (x P , Q 2 ) (4a) 



where 

A(A - 1) 8?r 



/ ( y ) = _ A ^- L ) ™ Re f d 2 h dz dz > e ik L (*'-*) 

y Jz>>z 



xp^(b,z,z') exp (4b) 

with ki = M y. The mean free path of the hadronic state X inside the nucleus is Lx = 
(&xn p(b, C)) > an d the assumption is made that the total cross section, &xn, for the 
state X with a nucleon is independent of the mass Mx - in which case we may take 
&xn ~ 20 — 30 mb, similar to the vector meson-nucleon cross sections. The uncertainty in 
the precise value of ctxn is not significant numerically for small and intermediate Q 2 , since 
there S^F 2 a is much smaller than S^ v 'F 2 a- At higher Q 2 , however, the larger nucleon sea 
gives a large P-exchange contribution, so that above Q 2 ~ 5 GeV 2 it starts to dominate. 

The combined effect of both the VMD and P-exchange mechanisms is a total shadowing 
correction that varies little for Q 2 between ~ 0.5 and 10 GeV 2 . In this transition region con- 
cerns have been expressed about possible double counting when adding the two components 
0. Since the VMD contribution is a higher twist effect, this problem could in principle be 
averted by keeping only the leading twist piece of F 2 p. Alternatively, the solution adopted 
by Kwiecinski & Badelek || was to exclude from the P-exchange process those final states 
X which have a mass comparable to that of the vector mesons {M x < M Xo ~ 1.5 GeV 2 ), 
and consequently take as the lower limit of integration in Eq.(f|a|) y m in = x(l + M X() /Q 2 )- 
Numerically, this will be important only at small x (x ^$ 10~ 3 ) and Q 2 , while for larger Q 2 
the separation becomes redundant because there y min — > x. 

In Figs. 1-3 we illustrate the effect of shadowing on the ratios of nuclear to deuteron 
structure functions (per nucleon) for various nuclei (with F 2 a = (Z / A)F 2p + (1 — Z/A)F 2n + 
S^F 2A + 5^F 2A ). ForF 2p and F 2n we use the parameterisations of Donnachie & Landshoff 



[ 15|] . (In Ref. |15| the neutron structure function was extracted from F 2 d without taking 
into account the possible shadowing corrections in the deuteron. However, these are much 
smaller than the shadowing corrections for the heavy nuclei considered here and would have 
little effect on the final A/D ratios.) 

Figure 1 shows the prediction of the model for the Xe/D structure function ratio, com- 
pared with the recent E665 data Q. The dashed curve is the result of the VMD model, 



while the solid curve includes in addition the P-exchange component. Clearly the bulk of 
the shadowing is due to the rescattering of vector mesons. Note that the data (and the 
curves) at each x correspond to a different average value of Q 2 , ranging from (Q 2 ) w 0.03 
GeV 2 at the lowest x value to (Q 2 ) ~ 10 GeV 2 at x ~ 0.1. The flattening of the curves as 
x — > is a direct consequence of the vanishing (like ~ Q 2 ) of the total F 2 a as Q 2 — > (i.e. 
a smooth extrapolation to the photoproduction limit). As expected, the role of two-particle 
correlations is found to be rather small, with the result of using an independent particle 
approximation differing only by some 5 — 10% from that in Fig.l. Of greater significance 
are the effects due to projectile attenuation, which reduce the amount of shadowing by up 
to 30 — 40%. Indeed, for larger nuclei the attenuation is even stronger, so that eventually 
saturation with A is reached. This is illustrated in Fig.2 where we plot F 2 a/ 1 F 2 d for v = 150 
and 250 GeV (at x = 0.0025), and compare with the recent Fermilab data for Ca,Xe and 



Pb, taken at 0.0018 < x < 0.0032 [16 



For completeness, in Fig. 3 we also compare the model predictions with the high-precision 
NMC data [g] on lighter nuclei (Ca and C), taken at 0.5 < Q 2 < 10 GeV 2 . Clearly the 
agreement with the data is quite good for both Ca and C. 

Having obtained good agreement with the data on heavy nuclei, we next calculate the 
shadowing corrections to the deuteron structure function in the region x ^ 10~ 5 and Q 2 <; 
0.03 GeV 2 , where the recent Fermilab data were collected 0. Traditionally in DIS on the 
deuteron, in which the proton and neutron are held together very weakly, nuclear effects 



have been ignored. Recent calculations [17,13 1 , however, have suggested small, but possibly 



non- negligible corrections at Q 2 ~ 4 GeV 2 . 

For 7*D scattering the Glauber formalism involves just the double scattering term in 
the multiple scattering series. In this case the shadowing corrections to F 2 d are given by 
Eqs.(|Ta|) and (Sal), except now the distribution of nucleons in the deuteron is given by the 



deuteron form factor P3| , Sofa ) = f£°dr(u (r) + w (r)) jo(|k|r), where u(r),w(r) are the 



S, -D-wave deuteron wavefunctions taken from realistic iViV potential models |Tj| . 



Apart from the VMD and P-exchange contributions, another potential source of shadow- 



ing in D is the exchange of (unreggeised) mesons, which was first shown by Kaptari et al. [19 
to give sizeable antishadowing corrections to F 2 d at Q 2 ~ 4 GeV 2 . In the non-relativistic 
approximation, this contribution is given by the convolution of an exchanged-meson dis- 
tribution fM/oiy) 0,0' w ith the meson structure function F2m(x/v,Q 2 )- Unfortunately, 
experimentally only F 2n is known, and then only in the region x > 0.2 and Q 2 > 4 GeV 2 . 
In addition there are no data on pion photoproduction which could constrain the Q 2 — > 
extrapolation parameters for F 2l ^. As a simple solution we take the (real) pion structure 
function measured in Drell-Yan reactions (as parameterised in Ref. | |20|| , with the fit that 



corresponds to a 20% pion sea at Q 2 = 5 GeV 2 ), and apply a phenomenological factor 
(Q 2 /(Q 2 + Ql)) , motivated by Regge theory. For Q 2 we use the same value as in the 
analysis of the ■yp data (see above), since Ql should be independent of the target hadron 
if it is to be interpreted as the 'mass' of the qq pair in the photon. Although the extrap- 
olation of F 2 m is the main uncertainty in the calculation, fortunately the meson-exchange 
contribution turns out to be small in the region of x and Q 2 covered by the E665 data 
(its contribution at x ~ 10~ 3 (10~ 2 ) is some 6% (13%) of the total shadowing correction to 
F 2 d)- This result is largely independent of the model D wavefunctions in Ref. [|T^] or the 
meson-nucleon form factor used. 

In Fig.4 we show the ratio of deuteron to proton structure functions with (solid curve) 
and without (dotted curve) the VMD, V- and meson-exchange shadowing corrections, and 
compare with the recent data from the E665 Collaboration ||. For each x the E665 data 
were taken at various photon energies, ranging from v ~ 300 GeV at the lowest x values 
down to v ~ 140 GeV at x « 10~ 2 . Also shown in Fig.4 are the NMC data at Q 2 = 4 
GeV 2 , and the corresponding theoretical curves (for x > 0.05). The net effect of shadowing 
is a decrease in F 2 d/F2 P by less than ~ 3% over the experimental range of x and Q 2 . As a 
fraction of the total F 2 d this amounts to < 2 — 3%. Consequently, for v ~ 150 — 350 GeV, 
the free neutron structure function is ~ 5 (2)% larger at x ~ 10~ 4 (10~ 2 ) than the bound 
neutron structure function, extracted from the deuteron data without account of shadowing. 
Unfortunately the E665 data at present do not distinguish between the shadowing and no- 
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shadowing scenarios, and it is essential that the errors be reduced before any unambiguous 
conclusions can be reached. As an independent check, it would be useful if the statistics 
and range of x of the neutrino-proton DIS data were also improved. This would allow 
the individual quark flavour distributions in the proton to be determined, and the neutron 
structure function inferred by charge symmetry. 

In conclusion, we have seen that the two-phase VMD/Pomeron-exchange model can give 
a quantitative description of the shadowing observed by the E665 and NM Collaborations 
in nuclear to deuterium structure function ratios at low x and Q 2 . The apparent saturation 
of shadowing at x ^ 10~ 3 in Xe/D is naturally reproduced by a careful extrapolation to 
the total structure function to the photoproduction limit. Furthermore, saturation with 
increasing A is found to be largely a consequence of the attenuation of the projectile debris 
inside the nucleus. The good agreement with the present nuclear data should also give us 
some confidence in the predictions for shadowing in the deuteron itself. 
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function parameterisations, and N.N.Nikolaev and A.Yu.Umnikov for helpful communica- 
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Figure Captions. 

1. Shadowing in the Xe nucleus. The dashed curve is the result of the VMD model, while 
the solid curve includes in addition the Pomeron-exchange contribution. The data are from 
the E665 Collaboration [Q], and both the statistical and systematic errors are shown. 

2. Shadowing as a function of mass number A. The theoretical curves are evaluated at 
x = 0.0025, for v = 150 GeV (dashed) and v = 250 GeV (solid), and compared with the 
preliminary E665 data for Ca, Xe and Pb, taken for 0.0018 < x < 0.0032 [lfj . 



3. Shadowing in the Ca and C nuclei. The data are from the NMC [[J, with statistical and 
systematic errors added in quadrature. 

4. Ratio of deuteron to proton structure functions with (solid curve) and without (dotted) 
shadowing corrections. The data (full circles) are from the E665 Collaboration (with statis- 
tical and systematic errors shown separately), and from the NMC (open circles) which were 
taken at Q 2 = 4 GeV 2 . 
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